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Abstract 

One and two loop self-energies are worked out explicitly for a heavy scalar field inter- 
acting weakly with a light self-interacting scalar field at finite temperature. The ring/ daisy 
diagrams and a set of necklace diagrams can be summed simultaneously. This simple 
model serves to illustrate the connection between multi-loop self-energy diagrams and 
multiple scattering in a medium. 
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The study of how and why the properties of particles change in a thermal medium 
is fascinating and of particular importance in the areas of cosmology, astophysics, and 
high energy nucleus-nucleus collisions. In this paper we investigate a very simple model 
involving a heavy scalar field with vacuum mass M interacting with a much lighter scalar 
field with vacuum mass m <^ M at temperature T -C M. The lighter field has a quartic 
self-interaction that is stronger than the interaction between the light and heavy fields. 
Although the model considered here is not realistic it has many properties in common 
with more physically interesting theories, such as the decay of a very heavy boson in the 
early universe or the emission of a very massive virtual photon in high energy nuclear 
collisions. The simplicity of this model allows a very clear mathematical evaluation of 
the heavy boson self-energy at finite temperature and its physical interpretation in terms 
of multiple scattering in the many-particle medium. One loop self-energies were studied 
and physically interpreted in a beautiful paper by Weldon [|I|. Some examples of two 
loop calculations are 0, 0, ^, |] but since those are oriented toward real physical appli- 
cations the resulting formulae are necessarily much more involved and oftentimes tedious 
to reproduce. We should also like to point out that at low to moderate densities the self- 
energy can always be expressed as an integral over the scattering amplitude for scattering 
from particles in the medium. Frequently one can use experimental input of scattering 
amplitudes to construct the self-energy 0. A formal discussion of this is was presented 
by Jeon and Ellis 0. 

The heavy scalar field is labeled H and the much lighter scalar field is labeled 0. The 
interaction Lagrangian is 

Ci = gHct)^ - \ct)\ (1) 

The cubic coupling g has the dimension of mass while the quartic coupling is dimensionless. 
It is assumed that <^ AM^ in order that we can ignore any self-energy diagrams 
involving an internal H line. This is analogous to the problem of the coupling of a photon 
to quarks where the electromagnetic coupling is much smaller than the strong coupling. 

The resulting one and two loop diagrams for the H are drawn in Fig. |l|. We are 
interested in the limit that the H is at rest with respect to the medium. (The interested 
reader may easily generalize the results that follow.) We will first evaluate the diagrams 
in Euclidean space and then analytically continue to Minkowski space. Denoting the Mat- 
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Figure 1: The one and two loop contributions to the H self-energy under the condition 
g^ ^ AM^. The thick dashed line is the if, the solid line is the 0. There is a factor g at 
each three point vertex and a factor —A at each four point vertex. 
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subara frequency of the external H line by Un = 27rnT we can find explicit mathematical 
expressions for these diagrams in the usual way [H]. 
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These summations may be evaluated directly or by the use of contour integration. 
The simplest summation is 



F(m,T) = T^ j 
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where uj = \/p'^ + m'^ and ?t.be(3;) = 1/(6"^ — 1) is the Bose-Einstein occupation number. 
The next simplest is 
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The final summation needed is 
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We can express the one and two loop self-energies in terms of these functions. 

Hi = -2g^G{u^,m,T) 
U2A = 24A/G2K, m,T) 

d 
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In Euclidean space the self-energy is real. Note that both the one and two loop contribu- 
tions can be expressed in terms of one dimensional integrals. 
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In order to obtain physical quantities, such as the dispersion relation, scattering, 
decay, and production rates, we must analytically continue to Minkowski space. The 
retarded self-energy is obtained by the replacement Un —iM + e with M > and 
e ^ 0+. Then 




(nBE(M/2r) + 1)^ - ni^{M/2T) (11) 



ImG = 

Here P.V. stands for principle value. The 1 following the Bose-Einstein occupancy number 
in the real part above represents the vacuum contribution. To really evaluate it we should 
go back to a four dimensional integration over momentum in Euclidean space and place an 
upper cutoff on the integration, or else analytically continue to 4 + e dimensions; in either 
case renormalization would proceed as usual. See the appendix. The term (?t,be + 1)^ in 
the imaginary part represents the decay H ^ (fxp with the Bose-Einstein enhancement 
factor in the final state; the term represents the production 00 H. 

The real part of 11 represents a shift in the mass squared of the H and the imaginary 
part represents the rate of its decay, scattering, or production. It is instructive to consider 
the hmit m -C T < M. 

F ^ h vacuum (14) 

12 ^ ^ 

The imaginary part of the H self-energy at the one and two loop order now takes a very 
simple form. 



im n = — 
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Stt 



(riBE + 1)' - 4e] {l + (4 - 2)A^ + ...| (15) 



This formula may be interpreted thusly. The overall sign results from the convention used 
for the self-energy. A minus sign indicates loss or decay of the H, a plus sign indicates 
gain or production. The factor (riBE + 1)^ shows the Bose-Einstein enhancement of the 
final state, the factor ril^ shows that the probability of production is proportional to the 
probability to find two of the requisite energy in the initial state. (The hbe are of 
course evaluated at the energy M/2.) The correction to order A has contributions from 
the diagrams 112^ giving the factor 4 and from 112^ giving the factor —2. Is there a simple 
interpretation of these results? 

The imaginary part can be obtained by cutting a diagram in half. Cutting the one 
loop diagram, as in Fig. ^ illustrates the physical process of decay of the H into 00. 
(For definiteness we focus on the absorption or decay processes. Reversing the in and out 
states gives the production rate.) 
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Figure 2: Cutting the one loop self-energy diagram. 



Cutting the diagram for 112^ does not yield the square of an amplitude. Instead it 
illustrates an interference between an amplitude of order zero and an amplitude of order 
A. See Fig. |^. Note that there should be a spectator with energy uj in the simpler half 
of the interference graphs as shown in previous work [1^]. We have not included it 
here to reduce the number of sub-figures within each figure that need to be drawn. The 
same also applies to Fig. ^ below. The loop that is not cut involves a Matsubara sum. 
This can be expressed as a vacuum loop plus a finite temperature piece. The former is 
just the order A correction to the vertex coupling the H to (f)(f). The latter can be viewed 
as the absorption of a from the heat bath with energy uj and momentum p and the 
emission of a into the heat bath with the same energy and momentum. The initial and 
final states are the same as in the amplitude of lower order in A provided a spectator 
is included in the simpler half of the interference graph as mentioned above, hence the 
reason that it can interfere. There are two diagrams with a single intermediate off- shell 
0. The kinematics are such that one of them has propagator l/q\ with = M{M + 2uj), 
and the other has propagator with = M{M — 2uj). Adding these two propagators 
gives 2/(M^ — 4cl'^). This is precisely the origin of the mathematical expressions displayed 
in eqs. (6) and (12) above. 
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Figure 3: Cutting the two loop diagram for Ti2A results in an interference between two 
amplitudes of different order in A. A loop involves a Matsubara sum unless otherwise 
indicated as a vacuum loop. The g+ and g_ are the four-momenta of intermediate states. 
Combinatoric factors are not displayed. 
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Figure 4: Cutting the two loop diagram for 1125 is just the second term in the Taylor 
series expansion for the one loop diagram shown in figure 2 with the bare propagator 
replaced by the dressed one. In effect the in and out (p lines should be dressed to represent 
collective excitations in the thermal medium rather than vacuum particles. 

Cutting the diagram for 1125 is shown in Fig. ^ Again there is an interference be- 
tween amplitudes of different order in A. One is tempted to interpret this as scattering 
of a particle from the heat bath with one of the outgoing from the H decay. However, 
because there is an intermediate line on-shell the diagram is not well-defined mathemati- 
cally (one should also read the paragraph near the end of the paper which discusses how 
in the imaginary-time formalism the contour integral allows one to avoid this problem 
altogether). This should be no surprise: one does not dress external lines when evalu- 
ating Feynman diagrams for scattering or decay. Indeed, this diagram really arises from 
an expansion of 11 in terms of a dressed propagator and bare vertices rather than an 
expansion in terms of a bare propagator and bare vertices. The one loop contribution to 
the self-energy of the is momentum and frequency independent and is simply given by 
12AF(m, T) . The dressed propagator is 1/ (a;| + + m^) where = rn? + 12AF(m, T) . 
One sees very clearly that I1.2B is the second term in a Taylor series, the first term be- 
ing Hi. All diagrams of this type can easily be summed by making the replacement 
ni(u;„, m,T) Iii{uJn,rnT,T). Physically what is happening is that the H is decaying 
into collective excitations with the quantum numbers of the rather than into vacuum- 
like excitations. In the context of gauge theories this goes by the name of hard thermal 
loops |jTl|. It is left as an exercise for the reader to show that the combinatoric factor 



for each diagram in this class is exactly reproduced by this selective resummation. Such 
perturbative diagrams go by the name of ring or daisy diagrams. One can go even further 
by summing super ring or super daisy diagrams by solving the integral equation 

rnl = m^ + 12XF{mT, T) (16) 

for the effective mass Examples of ring/daisy and super ring/daisy diagrams are 

shown in Fig. ^ (a) and (b) respectively. 
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Figure 5: (a) an example of a ring/ daisy diagram, and (b) an example of a super ring/ daisy 
diagram. 



There is another class of disgrams that are easily summable. We refer to these as 
necklace diagrams. They have exactly the same topological form as 112^ but with N loops 
instead of 2. The overall combinatoric factor is —12^/6. The expression for the one with 
N loops is 

12^ 

nnecklace ^ _ 2( _^)7V-1^A. ^^^^ y ) _ ^^7) 

O 

This is a geometric series readily summed to give 

_ 2g-'G{u^,m,T) 

J-J-necklace " 



l + 12XG{uJn,m,T)' ' ' 

Cutting the = 3 necklace, for example, leads to the diagrams shown in Fig. ^ and 
Fig. ^ At first one think that summing the necklaces amounts to a temperature coupling 
for H(f)(j), namely, g'^{T). This is not a good way to think because not only would this 
effective coupling have a dependence on the external Matsubara frequency u;„ but it will 
become complex when analytically continuing to Minkowski space. 



-ooo- 
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Figure 6: One way of cutting the = 3 necklace diagram. Combinatoric factors are not 
displayed. 
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Figure 7: Another way of cutting the = 3 necklace diagram. All the loops between the 
brackets are evaluated at T=0. Combinatoric factors are not displayed. 



One can further sum all diagrams of the ring and necklace types simultaneously simply 
by replacing m in the necklace formula with TtTt- 



n 



ring+necklace 



2g^G{un,mT,T) 
+ 12AG(cj„, mT,T) ■ 



(19) 



Neglecting the denominator, and expanding the numerator to zeroth and first order in 
A reproduces IIi and 1125, respectively. Keeping the numerator to zeroth order and 
expanding the denominator to first order reproduces T\.2A- 

We can also calculate the imaginary part for Fig. |l] from the scattering amplitudes 
shown in Figs. 2-4. This procedure was used in and described in more detail in |T^ . 
Rather than calculating the loop diagrams in Euclidean space, analytically continuing to 
Minkowski space, and then finding the imaginary part, this procedure uses the scattering 
amplitudes as input and does a thermal averaging over particles in the heat bath. In this 
procedure, every line carrying a distinct 4-momentum is integrated over phase space with 
the appropriate thermal weight. The integrations are of course constrained by a global 4- 
momentum conserving delta-function. The amplitudes are evaluated at tree level, except 
for the explicit vacuum loops. There are no finite temperature loops. The amplitudes 
are computed according to conventional vacuum field theory. In those cases where an 
outgoing hne is in effect dressed by the scattering of a particle from the heat bath, as in 
Fig. 4, then the amplitude will have to be dealt with more carefully; see below and also 
ref. |]TD[ for a general discussion. 
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Applying the scattering amplitude method to the current calculation of the decay of 
the heavy scalar field H at rest, one starts directly with the interference graphs on the 
second line in Fig. |^. A phase space integration measure is assigned to each independent 
incoming and outgoing line except, of course, the H line. The two lines carrying the energy 
io will have only one and not two integral measures because the same 4-momentum flows 
through them. They come from the heat bath and go back into it with the same momen- 
tum, therefore they are not constrained by the overall energy-momentum conservation. 
The contribution from Fig. |^ to the decay rate F can be written as 



2i?F(Fig, 3) 



X 



l+r2BE(rf/T))(l+nBE(p^/T))(A^L,®A^(Fig.H)+h. c.) . (20) 



The factor of minus one-half is a symmetry factor. Here P = {E, P) is the 4-momentum 
of the H. In this paper, though, we specialize to the H at rest so that P = and P° = M. 
For the simpler half of the interference graph, Ai is just the tree amplitude 

Mtree = 1g . (21) 

The other half is more complicated. 

M,r.,:, = 2!x4!xj(-A)/p^(l + nB.(./r))(-^-^) 

= + {T = pai't} + 0(1/Af') . (22) 

The factors 2! and 4! come from the Hcf)^ vertex and the (f)^ vertex, respectively. The 
(1/2-1- nBE) is the distribution assigned to boson lines which correspond to absorption 
and emission of the same 4-momenta within the same Feynman graph according to one 



of the rules given in |T0[. Its origin is really just the same as Eq. (5). The denominators 
and introduced earlier come out automatically from calculating the last two graphs 
in Fig. 1^. They have a minus sign because of our conventions, namely, the propagator 
is positive in Euclidean space [§]. With Eq. (|2T|) and Eq. (p2D , Eq. ( pO|) can easily be 
evaluated to give 



2Mr(Fig, 3) = — + {T = Q part} -|- {non-leading terms} . (23) 



The other contribution to the H decay comes from the processes shown in Fig. ^. 
This time there is an intermediate line which is on-shell because it carries the same 4- 
momentum as an outgoing line. For that line, with 4-momentum pi, the phase space 
integral measure must be written out as a 4-momentum integral measure because there is 
a differentiation of the mass-shell constraining delta-function. The intermediate on-shell 
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line is mathematically absorbed into this differentiation of the mass-shell constraint and 
does not appear anywhere else in the calculation. In terms of the self-energy loop diagrams 
shown in Fig. 1, this derivative of a Dirac delta-function has its origins in a double pole 
arising from a self-energy insertion. The contribution from Fig. ^ can be written as 

2Mr„, , . 1 X 2 / / 2..'«(P?) (2^rS">(P - P. - P.) 



X 



1 + nBE(rf/T)) (1 + nBEipVT)](Mt,, ® A<(Fig. 0) + h. c.) . (24) 



The factor of two here is for the emission and absorption to occur on either of the two 
outgoing lines. A^tree is the same as above in Eq. (^I]) . The symbol S means to 
differentiate the Dirac delta-function with respect to its argument, taking the energy 
positive. The remaining half of the amplitude which contains the emission and absorption 
is 

r (Pk /I 

-M(Fig. 4) = 2! X 4! X g{-\) J j^-^ (- + nBE(^/r) 

= -2^AT2 + {T = part} . (25) 

The factors of 2! and 4! have the same origin as before, and the emission/absorption 
is associated with a distribution of (1/2 + ube)- The integration in Eq. ( p^ has to 
be evaluated with some care in order to manipulate correctly the differentiation of the 
mass-shell constraint. The result is 



2Mr(Fig. 4) = — ^ + {T = part} + {non-leading terms} . (26) 



So we have ^ ^ 

2M(r(pig. 3) + r(pig. ^ (2 - 1) + {T = part} . (27) 

The result from Eq. ([15|) can be recovered with 

r = r(Fig. 2) + r(Fig. n + r(Fig. ^ (28) 

using the identification 

2Mr (exp(-M/T) - l) = 2 Jm H . (29) 



One point in the calculation that should be elaborated concerns self-energy insertions. 
In terms of scattering amplitudes it amounts to dressing an external line, thereby causing 
a singularity on an internal line. The usual Feynman rules for scattering theory say that 
external lines should not be dressed but should be taken into account in the phase space. 
(In our example, this is the intermediate line in Fig. ^) This usually leads to the need 
to include at least part of the thermal self-energy in the field, and therefore changes the 
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mass-shell of the particle. This was shown quite elegantly by Weldon |T2[ for fermions. 
Donoghue and Holstein have done the same but in a Lorentz non-covariant manner |T3 . 
When calculating loops in the imaginary-time formalism, there is a rigorous solution 
to this problem. It is well-known that the contour integral method has a well-defined 
way of dealing with multi-pole singularities. The residue of the multi-pole is found by 
differentiating the non-pole part of the integrand by one power less than that of the order 
of the multi-pole. So, if one insists on no resummation of any sort, including those of 
super ring/daisy diagrams mentioned above, or even partial resummation, such as those 
of Weldon, Donoghue and Holstein, then the contour method is a good solution to the 
above problem. This method was also used in |10 



In conclusion we have accomplished several goals. First, we have explicitly evaluated 
the two loop self-energy for a massive boson weakly interacting with lighter bosons in 
thermal equilibrium among themselves. This was done first in the imaginary time for- 
malism of Matsubara. There is no ambiguity in analytically continuing to real time or 
real energies, as first proven by Baym and Mermin 40 years ago Second, a physical 



interpretation was found for all diagrams considered. Some of these involve Bose-Einstein 
enhancement of strongly interacting particles in the final state of the heavy boson decay. 
Some of them involve absorption and emission of particles from the heat bath with the 
same energy and momentum. Such processes may interfere with processes of other order 
in the stronger coupling A. Third, it was seen how important it is to expand in terms 
of resummed or interacting propagators. This indicates the importance of identifying the 
collective excitations of the system. These are all features that must be understood and 
taken into account when considering the kinetic theory of quark-gluon plasma in heavy 
ion collisions or in the early universe, for example. Then the role of quarks and gluons 
are analogous to the cf) bosons, and the role of the H is analogous to the electromagnetic 
field. 
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Appendix 



In this appendix we compute the vacuum self-energy of the H boson to one and two loop 
orders. The diagrams are those shown in Fig. 1. We use dimensional regularization (see, 
for example, chapter 4 of |]T5|) and take the mass m of the boson to be zero for simplicity 
of presentation. The two integrals which are needed are the F and G as displayed in eqs. 
(2-10). In Euclidean space: 



and 



2\2-uj 



(27r)2^ P 



1 



(30) 



(31) 



(27r)2^ P (l-Py 

Here / is the loop momentum and P is the external momentum, both Euclidean. The 
number of dimensions is 2a; with 2 — u; = e ^ O"*", and fi is the dimensional parameter 
which is arbitrary. Dimensional regularization gives F^ac = 0. 
Using the Feynman parametrization one arrives at 



+ 



+ 



(47r)2 
4 



p2 




^-2^E + \ll + (2 - 1e) In 




(32) 



-P2 and 



Analytic continuation to Minkowski space involves the replacements - 
In {4iT^'^/P'^) In {An fi"^ / P"^) + in. This results in Gvac having both real and imaginary 
parts. 

The imaginary part of the H vacuum self-energy can easily be found from the above 
results in conjunction with eqs. (8-10). At one loop order 



9' 



in 



and at two loop order 




2-fE + 2 In 




(33) 



(34) 



This needs a counterterm added to the Lagrangian to render it finite in the limit e — 0. 
We have the freedom to choose the counterterm such that the order A correction vanishes 
on the mass-shell. If we make this choice then 




(35) 
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where M is the physical mass of the H boson. In effect this defines the normahzation 
point of the couphng g of the H to two mesons. The sum of one and two loops is 



Im IIv 



57r 



3A, (M^' 



8^1 + filn(^) Stt 



(36) 



This involves the renormalization group running coupling. It displays infrared freedom, 
2^p2) ^ as P2 _^ 0, and the Landau pole at = M2exp(7rV3A). 

The above results can also be obtained from the vacuum diagrams shown in Figs. 2 
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to 4. The methods for doing this are well known: see, for example, chapter 2 of [^]. This 
is left as an exercise for the reader. 
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